Eigenfunctions of electrons in weakly disordered quantum dots: 
Crossover between orthogonal and unitary symmetries 
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A one-parameter random matrix model is proposed for describing the statistics of the local 
amplitudes and phases of electron eigenfunctions in a mesoscopic quantum dot in an arbitrary 
magnetic field. Comparison of the statistics obtained with recent results derived from first principles 
within the framework of supersymmetry technique allows to identify a transition parameter with real 
microscopic characteristics of the problem. The random-matrix model is applied to the statistics of 
the height of the resonance conductance of a quantum dot in the regime of the crossover between 
orthogonal and unitary symmetry classes. 
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I. INTRODUCTION 

Statistical properties of the electron eigenfunctions in 
disordered quantum dots have recently become a sub- 
ject attracting considerable theoretical and experimental 
interest - (|]. One of the reasons is that the prob- 
lem of particle motion in a bounded disordered poten- 
tial comprises a particular case of general chaotic sys- 
tems, such as quantum billiards 0. On the other hand, 
the development of the "microwave-cavity" technique 
- H has opened up a unique possibility for experimen- 
tal studies of wave-function statistics. Furthermore, the 
mesoscopic conductance fluctuations in a quantum dot in 
the Coulomb blockade regime contain fingerprints of the 
statistics of the eigenfunctions of confined electrons M, 
JlO| - [|l3| . This also provides experimental access [ fj"4|]l5| 
to the wave-function statistics. 

In order to characterize the eigenfunction statistics, we 
introduce two distribution functions. The first refers to 
the local amplitudes of electron eigenfunctions ip n inside 
a dot of volume V, 



P(v) = (s (v-V\^ n (r)\ 2 )) 
while the second distribution function 
Q (s) = {6(s- s'mifn (r))) 



(1) 



(2) 



is associated with the eigenfunction phases, <p n = arg-0„. 
Angular brackets (...) denote averaging over bulk impu- 
rities inside a quantum dot. 

When all electronic states are extended (metallic 
regime), the random matrix theory (RMT) |L6| predicts 
the Porter-Thomas distribution of eigenfunction inten- 
sity that only depends on the fundamental symmetry of 

The parameter (3=1 for a system having time-reversal 
symmetry (orthogonal ensemble), whereas (3 = 2 when 



time-reversal symmetry is broken (unitary ensemble), 
and (3 = 4 for a system having time-reversal symmetry 
and strong spin-orbit interactions (symplectic ensemble) . 
The RMT predictions are known to be valid in the limit of 
large conductance g = E c /A 3> 1 (here, E c = hD/L 2 is 
the Thouless energy, D is the classical diffusion constant, 
L is the system size, and A is the mean level spacing) 
providing the system size L is much larger than the elec- 
tron mean free path I. Corrections to the Porter-Thomas 
distribution calculated in the framework of the c-model 
formalism pj] are of order g^ 1 in the weak- localization 
domain. 

The limiting case of orthogonal symmetry corresponds 
to systems with pure potential electron-impurity scatter- 
ing, while the case of unitary symmetry applies to sys- 
tems in a strong magnetic field that breaks the time- 
reversal symmetry completely. For intermediate mag- 
netic fields, a crossover occurs between pure orthogonal 
and pure unitary symmetry classes. This crossover in 
the distribution function V (v) was previously studied 
within the framework of supersymmetry techniques j^j 
and within an approach fij exploiting the analogy with 
the statistics of radiation in the regime of the crossover 
between ballistic and diffusive transport. 

We consider below the problem of eigenfunction statis- 
tics of chaotic electrons in a quantum dot in an arbitrary 
magnetic field in the framework of RMT. Our treat- 
ment is related to the case of the metallic regime, and 
describes statistical properties of eigenfunctions ampli- 
tudes and phases in the regime of the orthogonal-unitary 
crossover. The results are applied to describe distribu- 
tions of level widths and conductance peaks for weakly 
disordered quantum dots in the Coulomb blockade regime 
in the presence of an arbitrary magnetic field. 



the quantum dot: V (v) = {^j ' exp ( - <r/2 ). 
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II. DISTRIBUTION OF LOCAL AMPLITUDES 
AND PHASES OF EIGENFUNCTIONS 



In order to study the statistical properties of electron 
eigenf unctions within the RMT approach, we replace 
the microscopic Hamiltonian 7i of an electron confined 
in a dot by the N x N random matrix H = SpeSZ 1 
that exactly reproduces the energy levels e n of the elec- 
tron in a dot for a given impurity configuration. Here, 
e =diag(ei, ...,£jv), Sp is a matrix that diagonalizes ma- 
trix H (parameter j3 reflects the system symmetry), and 
it is implied that N — > oo. An ensemble of such random 
matrices reproduces the electron eigenlevels for micro- 
scopically different but macroscopically identical realiza- 
tions of the random potential, and therefore it should 
describe the level statistics of electrons inside the dot. 
As long as we consider the metallic regime, the relevant 
ensemble of random matrices is known to belong to the 
Wigner-Dyson class ■ Correspondingly, averaging 

over randomness of the system is replaced by averaging 
over the distribution function P [H] cx exp {— trV [H]} of 
the random-matrix elements, where V (e) is a so-called 
"confinement potential" that grows at least as fast as 
|e| at infinity pqjl9|| . In such a treatment, the columns 
of a diagonalizing matrix Sp contain eigenvectors of the 

1/2 

matrix H, that is ipj (i - ;) = (N/V) (Sp)^, provided 
the space inside a dot is divided onto N boxes with 
radius vectors r$ enumerated from 1 to N (the coeffi- 

1/2 

cient (N/V) is fixed by the normalization condition 
/ y dr|^-(r)| 2 = l). 

For the case of pure orthogonal symmetry, the eigen- 
functions are real, being the columns of the orthogonal 

1/2 

matrix Si (up to a normalization constant (N/V) ). 
When time-reversal symmetry is completely broken, 
the eigenfunctions are complex and may be thought 
of as elements of the unitary matrix 5*2 . We note 
that in this case, the real and imaginary parts of 
the eigenfunction are statistically independent, and 

(\~Reipj (r)| 2 ^) = /|lm-0j (r) | 2 ^> - It is natural to as- 
sume that in the transition region between orthogonal 
and unitary symmetry classes, the eigenfunctions of elec- 
trons can be constructed as a sum of two independent 
(real and imaginary) parts with different weights, such 



as |Im^-(r)r / |Re^(r 



= 7 



where the transi- 



tion parameter 7 accounts for the strength of the sym- 
metry breaking. This means that in the crossover region, 

1/2 

the eigenfunctions (V/N) 1 ipj (r^) can be imagined as 
columns of the matrix S, 



5 = 7TT^(° + Z7 °)' °° T = 1 ' OO t = 1, 



(3) 



composed of two independent orthogonal matrices O and 
O. The parameter 7 in the parametrization Eq. (j^) 
governs the crossover between pure orthogonal symme- 
try (7 = 0) and pure unitary symmetry (7 = 1). The 
values < 7 < 1 correspond to the transition region 
between the two symmetry classes. From a microscopic 
point of view, this parameter is connected to an external 
magnetic field. (This issue will be discussed below.) 
From Eq. (0), we obtain the moments 



M7,A0 = (|£ y -| 2p ) 



1 



s (I + 7T 



q=0 



O; 



2(p-9) 



(4) 



Here (. ..) n stands for integration over the orthogonal 
group pOj . The summation in Eq. (^) can be carried 
out using the large- N formula 



leading to 

M P (7, N) ■■ 



l 2 P\ 



2_y T( P 

N J 



p^PAX) x= i, 

NP XP ' 2 17 7 



(5) 



(0) 



which is valid over the whole transition region < 7 < 1 
in the thermodynamic limit N — > 00. (Here P p (X) de- 
notes the Legendre polynomial). One can easily see from 
the properties of Legendre polynomials that the ansatz 
introduced by Eq. (||) leads to the correct moments in 
the both limiting cases 7 = \p. p (0, N) = (2p - l)U/N p ] 
and 7= 1 [fi p (l,N)=p\/NP]. 

Using the RMT mapping described above, we can write 
the distribution function V (v), Eq. ([!]), in the form 



V(v) 



exp ( 



-iujv 



(iujNf 



-'^r(p + i) 

whence we get, with the help of Eq. @ 



l 2 P\ 



, (7) 



(8) 



where Io is the modified Bessel function. Equation (||) 
gives the distribution function of the local amplitudes of 
electron eigenfunctions in a quantum dot [see Fig. 1]. It 
is easy to confirm that for pure orthogonal (X — > 00) and 
pure unitary (X — 1) symmetries, Eq. (ra) yields 



V (v) = Xexp (-vX 2 ) I (vX\/X 2 -l) , 



V(v)\ 



1 



x- 



V(v)\ 



/2^ 6XP V 2 



exp (— v) . 



(9) 
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A different approach to the issue of eigenfunction 
statistics was recently proposed in Ref. |Q, whose au- 
thors obtained a single-integral representations for V (v) 
[see their Eq. (17)]. It can be shown that our Eq. (||) 
coincides with Eq. (17) of the cited work provided the 
parameter A is related to the parameter Y appearing in 
Ref. § by A = y/T+T/Y. 

The advantage of the random-matrix approach is that 
it allows one to calculate in a rather simple way, the 
distribution functions for many other quantities, in par- 
ticular, the phase distribution given by Eq. (|^). 

Since within the framework of the proposed random- 
matrix model sin ip j (r) 



7<V 



0% 



7 2 2 , and 



(sin 2p+1 ipj (r))_ ~ = 0, we have to compute 



Q(s) 



du> 
2^ 



exp(- 



-iws) 

p=0 



(-1) P (7^) 

r(2p + i) 



2p 



-o-2p(7). (10) 



where the average 
P2p (7) = 



O 



2p 



H+7 2 ^.) 



(11) 



0,0 



can be calculated using its integral representation, 



CT 2p (7) 



r(p) 



_1 exp |— x C 



x(6ff I dxxP- 1 — * ~ 1 " 2 ' ~ 2 " 2 



)} 



0,0 



(12) 



Expanding the exponent in the integrand of Eq. ( |l2| ) 
yields 

_ » (-!)« a* 



a2pi7) = W)L dxx Sn^TT 



E 

k=0 



,2k 



Oi 



2{p+q-k) 



(13) 



Using Eq. (§), we obtain after straightforward calcula- 
tions 



°2p (7) = 



1 



rfA 



ir Jo V\(i + \)h 2 + xf 



(14) 



Equations (|1C|) and (|14| ) yield the following formula for 
the phase distribution function in the crossover regime: 



7 



1 



^2 7 2 + S 2 (l- 7 2 )- 



(15) 



As can be seen from Eq. (|l_5j), the limiting case of pure 
orthogonal symmetry is characterized by the 5-functional 
phase distribution 



CM 



7^0 



(16) 



whereas the case of pure unitary symmetry is described 
by the uniform distribution 



Q(s) 



1 

2tt' 



(17) 



In the crossover region, the calculated phase distribution 
1 ( s ) — |vl — s 2 Q (s) displays a smooth transition from 
(5-functional distribution in the case of pure orthogonal 
symmetry (X — > 00) to the uniform distribution in the 
case of pure unitary symmetry (A = 1), see Fig. 2. 

In order to relate the phenomenological parameter 
7 (or A) entering the parametrization of Eq. (^) to 
the magnetic field breaking time-reversal symmetry in 
the real microscopic problem, we compare the distribu- 
tion given by Eq. (|8|) with the exact distribution de- 
rived within the framework of microscopic supersymme- 
try model (Eq. (13) in Ref. Although these distri- 
butions have different analytical forms, good numerical 
agreement is observed even for the tails of distribution 
functions after appropriate rescaling of our phenomeno- 
logical parameter X [see inset in Fig. 1, where the dis- 
tribution function ip (r) = ItV (t 2 ) for different values 
of parameter A is plotted] . Analysis of the behavior of 
(p (t) in both cases in the region of small r yields the 
following relation between transition parameter A and 
microscopic parameter X m appearing in Ref. 0: 



A = 2Xi 



exp 1 



(X 2 m ) \ $! (A m ) [Ai (A m ) - A 2 (X m )] + A 2 (A m ) - 



1 - $1 (X n 



-Ai (A„ 



(18) 



Here $ x (A m ) = [exp {-A 2 n } /X m ] f* m dyexp {y 2 }, and 

A n (A m ) 



on+l v2ri+l 



(19) 
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In accordance with Ref. |2j, the microscopic parameter 
X m = \4>/<f>o\ (a g E c / A) 1 ^ 2 , where a g is a factor depend- 
ing on the sample geometry, cj> is the magnetic field flux 
penetrating into the cross-sectional area of a sample, 
and cj>o is the flux quantum. It can be seen from Eq. 
( |l8| ) that for very weak breaking of time-reversal sym- 
metry (X m <C 1), X w 2y/ir/3X m . In the opposite limit 
of weak deviation from unitary symmetry [X m 3> 1), 
X* 1 + 1/2X1, 

The connection between parameters X and X m given 
by Eqs. ( |l8| ) and (|l9|) allows to use the simple random- 
matrix model proposed here for describing various phe- 
nomena occurring in quantum dots in arbitrary magnetic 
fields. 



III. STATISTICS OF RESONANCE 
CONDUCTANCE OF A QUANTUM DOT 

The issue of eigenfunction statistics is closely con- 
nected to the problem of conductance of a quantum dot 
that is weakly coupled to external leads in the Coulomb 
blockade regime. At low temperatures, the conductance 
of a dot exhibits sharp peaks as a function of the ex- 
ternal gate voltage. The height of conductance peaks 
strongly fluctuate since the coupling to the leads depends 
on the fluctuating magnitudes of electron eigenfunctions 
near the leads. Thus far both theoretical jlO| - []l3| and 
experimental 0,^5| studies were restricted to the pure 
symmetry classes (with conserved or completely broken 
time-reversal symmetry) even in the simplest case of two 
pointlike leads. Here we present the analytical treatment 
of the problem for arbitrary magnetic fields. 

Let us consider non-interacting electrons confined in 
a quantum dot of volume V with weak volume disorder 
(metallic regime) . The system is probed by two pointlikc 
leads weakly coupled to the dot at the points r L and r R . 
The Hamiltonian of the problem can be written in the 
form Q 

W ^(' :V + S A ) 2 + C, M 



TH 



-V [a R S (r - r R ) + a L 5 (r - r L )] 



(20) 



where U consists of the confinement potential and the po- 
tential responsible for electron scattering by impurities, 
th is the Heisenberg time, and cxr^l) is the dimensionless 
coupling parameter of the right (left) lead. We also sup- 
pose that the coupling between the leads and the dot is 
extremely weak, a R (i,\ <C 1, so that the only mechanism 
for electron transmission through the dot is tunneling. 

It can be shown that the heights of conductance peaks 
are entirely determined by the partial level widths 



2 a.R(£) T , 1 , / \|2 

1vR(L) = —V | Vv {r R {L) ) I 

7~H 



(21) 



in two temperature regimes. When T <C o:.r(.l)A, the 
heights g v — +%G V are given by the Breit-Wigner for- 
mula 



9v 



{luR + Iul) 



(22) 



At higher temperatures, a R (nA « T « A, the heights 
are determined by the Hauser-Feshbach formula 



9u 



2T 7„_r + 7„l 



(23) 



Using the distributions Px (ivRIL)) 01 the partial level 
widths, Eqs. (H|) and Eq. (|), 



p x \1uR(l) ) = — exp 



2a 



R{L) 



X 2 t h 



2a 



R(L) 



-JvR(L) 



T ( xVW^Ith 

X J o luR{L) 



(24) 



and assuming that the eigenfunctions of electrons near 
the left and right contacts are uncorrelated, that is 
\r R — 3> A, we can derive analytical expressions for 
the distribution of the conductance peaks heights. 



A. Breit-Wigner regime 

From Eqs. (p^), we conclude that the distribution of 
the conductance peaks is given by 



/>oo />oo 

Rx {9u) = / dlvR j AluhPx {Ivr) Px (lui 
Jo Jo 



x5[g v - 



^IvRlvI 



{ivR + IvlY 



) 



(25) 



Straightforward calculation lead to the following integral 
representation 



Rx {g v ) = 



e (1 - g v ) 
V 1 - 9u Jo 




fc=(±) ' J 

(26) 
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where 



and with conserved time-reversal symmetry: 



= 1 ± yi^;, 

/ (±) = a± Va 2 - 1, 



a= 2 





(27) 



The integral in Eq. (£26|) can be calculated, yielding 
the distribution function for the conductance peaks 



Rx {g v ) 



E 



X e(l-g„) 
2ir VI - g v 



where the following notation has been used: 

Mi = a + ^/T^gZVa 2 - 1, 



(28) 



fcj = v 7 ^ 2 - 1 



(29) 



and E (fc) stands for the elliptic integral. Equation (EE 
is valid for arbitrary magnetic field flux, and it interpo- 
lates between the two distributions corresponding to the 
quantum dot with completely broken time-reversal sym- 
metry 

^..W- t'i^'t"-*^ (30) 



2^1^; [l + g v (a?-ij\* 
T*(£)= £exp(-aX 2 £)£ 



Rx^oo {gv) — 



e(i-g„) a 

^VgW 1 - ~gl 1 + 9v (° 2 - 1) 



(31) 



Comparing Eqs. ( ]2q ) and (|31|), one can see that the 
influence of the magnetic field is most drastic in the re- 
gion of small heights of the conductance peaks, g v ~ 0, 
[Fig. 3]. Equations (§§) yields 



Rx(0) 



X 



2a' 



1 



(32) 



Thus, we conclude that the probability density of zero- 
height conductance peaks decreases from infinity to the 
value (2a 2 — l) ypH jZX m at arbitrary small magnetic 
fluxes penetrating into the sample. 



B. Hauser-Feshbach regime 

When a#(£)A <C T <C A, the distribution of the con- 
ductance peaks is given by 

/>oo />oo 

F x {g v ) = / d<y uR / d^ vL P x (j vR ) P x (riuh) 
Jo Jo 



x5[g v 



7T luRlvL 
2T IvR + 7kL 



(33) 



Carrying out the double integration and using Eq. (|24|), 
we obtain the following formula for the function T x (£) = 
Fx {gv) Tfy/otijOiR/ 2Tt jj , where the dimensionless vari- 
able £ = 2g v TTH I TT^otLUR: 



dfj, 1 + 



exp 



(34) 



This equation interpolates between two limiting distribu- 
tions 



7x(£) = 4r£exp(-aX 2 £) 



T x= i (0 = £ exp (-«£) [tf (£) + aK, (£)] (35) 



and 



for unitary and orthogonal symmetries, respectively. 

The distribution function Eq. (|m|) can be represented 
as a series 



oo oo 



yy - 



n— m— 



XV^ 2 - 1 



2(n+m) 



f(+) 



2(n— m) 



s-l+2m / _ , \ \ 

(2{n + m + l)\ I (+) 



s=-l-2r; 
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that is convenient for the calculation of the distribution 
of conductance peaks in the case of weak deviations from 
unitary symmetry. As is the case of very low tempera- 
tures, the distribution function of the heights of the con- 
ductance peaks given by Eq. (|34|) is most affected by 
magnetic field in the region g v ss 0: 

T x (0) = aX, (38) 

and therefore, the probability density of the conductance 
peaks of zero height immediately drops from infinity to 
the value Tx (0) ~ 2a^/n/3X m when an arbitrary small 
magnetic field (X m <g; 1) is applied to a system. 

IV. CONCLUSIONS 

We have introduced a one-parameter random ma- 
trix model for describing the eigenfunction statistics of 
chaotic electrons in a weakly disordered quantum dot 
in the crossover regime between orthogonal and unitary 
symmetry classes. Our treatment applies equally to the 
statistics of local amplitudes and local phases of electron 
eigenf unctions inside a dot in the presence of an arbitrary 
magnetic field. The transition parameter X entering our 
model is related to the microscopic parameters of the real 
physical problem, and therefore the distributions calcu- 
lated within the proposed random-matrix formalism can 
be used for the interpretation of experiments. 

This random matrix model has also been applied to 
describe the distribution function of the heights of con- 
ductance peaks for a quantum dot weakly coupled to ex- 
ternal pointlike leads in the regime of Coulomb blockade 
in the case of crossover between orthogonal and unitary 
symmetry. We have shown that the magnetic field ex- 
erts a very significant influence on the distribution of the 
heights of the conductance peaks in the region of small 
heights. The effect of the magnetic field consists of re- 
ducing the probability density of zero- height conductance 
peaks from infinity to a finite value for arbitrary small 
magnetic flux. 
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Fig. 1. Distribution function (p (r) for different parameters X. Inset: the same distribution function calculated 
from Eq. (||) (lines) and from Eq. (13) in Ref. (figures) for parameters X and X m connected by Eq. ([l8|). Solid 
line: X = 2.56, triangles X m = 0.5; dashed line: X = 1.52, squares: X m = 1; dot-dashed line: X = 1.14, circles: 
X m = 2. 

Fig. 2. Phase distribution q (s). Parameters: X = 1.5 (curve 1), X = 5 (curve 2), X = 10 (curve 3). 
Fig. 3. Distribution of the heights of the resonance conductance peaks in the Breit-Wigner regime, a = 1. 
Parameters: X — 1.2 (curve 1), X = 2 (curve 2), X = 4 (curve 3). 
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